function y=x_m1(x);

dim=size(x);
n=dim(1);

%y=[ones(n,1) x(:,1)];

%y=[ones(n,1) x];

%y=[ones(n,1) x x.^2 x(:,1).*x(:,2) x(:,1).*x(:,3) x(:,2).*x(:,3)];

y=[ones(n,1) x(:,1)+x(:,2)+x(:,3)];